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Abstract 

The staircase model is a recently discovered || one-parameter family of 
integrable two-dimensional continuum field theories. We analyze the novel 
critical behavior of this model, seen as a perturbation of a minimal conformal 
theory M p : the leading thermodynamic singularities are simultaneously gov- 
erned by all fixed points M p , M p _i, . . . , M3. The exponents of the magnetic 
susceptibility and the specific heat are obtained exactly. Various corrections 
to scaling are discussed, among them a new type specific to crossover phe- 
nomena between critical fixed points. 
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In the past few years, there has been renewed interest in exploring the conse- 
quences of integrability in two-dimensional statistical systems. At a critical point, 
an infinite number of integrals of motion appears if the theory is conformally in- 
variant; one has even a partial classification of such universality classes ]l|]. For 
some perturbations away from criticality, a subset of these integrals survives and 
makes the theory solvable even at finite correlation length £. Two types of such 
systems with a generic (p — l)critical point (p = 3, 4, . . .) described by the minimal 
conformal theory M p have been known: 

(a) lattice models 0, whose manifold of integrability is parametrized about the 
fixed point M p by a relevant temperature-like thermodynamic parameter as well as 
marginal and irrelevant parameters governing the lattice effects, and 

(b) exact factorizable scattering matrices for certain relevant perturbations of the 
conformal fixed points M p |J. These define massive continuum field theories that 
describe universal scaling behavior. The thermodynamic Bethe ansatz is a way 
to calculate the universal ground state energy of the associated Hamiltonian on a 
circle of circumference R, 

E (R,O = ^f(R/O- (i) 

The scaling function f(p) shows a simple crossover from the thermodynamic regime 
R 3> £ to the conformal regime R <C £; its ultraviolet limit is determined by the 
central charge c of the asymptotic conformal theory: lim p _>.o f(p) — — c/12. 

The "staircase model" is a one-parameter family M(9q) (9 > 0) of factoriz- 
able scattering theories with a new and more intricate scaling behavior, discov- 
ered very recently by ALB. Zamolodchikov ||. These theories contain a single 
type of massive particles that are characterized by the purely elastic S-matrix 
S(9, 9q) = (sinh 9 — i cosh 26 Q ) / (sinh 9 + i cosh 2# ), written in terms of the Lorentz- 
invariant rapidity difference 9. The scaling function /(p, O ) shows a staircase pat- 
tern that interpolates between all central charges c p (see fig. 1). Hence the theory 
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M(9q) is described by a renormalization group trajectory that comes close to each 
fixed point M p for a RG "time" interval 8q, whereafter it crosses over to the next 
lower fixed point M p _i [|J. 

In another paper |J , I argued that from a Lagrangian point of view, this model 
can be understood as any fixed point theory M p , perturbed by a linear combination 
of its weakest relevant scaling field </>(i,3) and its leading irrelevant scaling field 0(3, i), 

£ = £p + t p 0(l,3) - £p0(3,l) > (2) 

t p > and i p > are dimensionful coupling constants with renormalization group 
eigenvalues y p — 4/(p+ 1) and f/ p = — 4/p, respectively The purpose of this Letter 
is to show in detail that this leads to a unique critical behavior as a function of 
the relevant "temperature" t p : for t p = 0, the theory is in the universality class of 
M p , but the thermodynamic singularities as t p ~\0 are determined not by M p alone, 
but simultaneously by all fixed points M p , M p _i, . . . , M 3 . Specifically, I obtain the 
exact exponents of the magnetic susceptibility 

x(tp) ~ t; 7p (3) 

and the specific heat 

C(tp)~t; a »iogtp (4) 

and discuss various corrections to scaling. 

Consider first a crossover between two critical fixed points which is characterized 
by a renormalization group trajectory u(r). The two-point correlation functions 
Gij(r,u) = ((f>i(Q,u)(f)j(r,u)) of local fields (pi are path-ordered integrals 

G ii (e* a - 7 V,u(7i)) = 

G kl (r,u(r 2 )) (PeM- jTdrx^Cr)))^ (Pexp(- jT'drx^^)))) (5) 

over the matrix of anomalous dimensions x(w), which is given in perturbation theory 
about the ultraviolet fixed point by 

xf(u)=x i i (0) + C i ' h u k + O{u 2 ) (6) 
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in terms of the asymptotic anomalous dimensions and operator product coefficients 
C? k . These equations imply the mixing of fields under the renormalization group: 
the infrared asymptotic scaling fields, i.e. the eigenvectors of the matrix x(u*), are 
linear combinations of the ultraviolet scaling fields, i.e. the eigenvectors of x(0). 
If e.g. the local magnetization o~(r) is expanded in the infrared scaling basis, the 
lowest appearing eigenvalue x gives its scaling dimension and the other eigenvalues 
x', . . . contribute corrections to scaling. Even if a(r) is chosen a pure scaling field 
in the ultraviolet, the mixing of fields dictates that the infrared asymptotics of its 
two-point function G a (r) is G a (r) ~ r ~ 2:c (i _|_ Ar~ 2 ( x '~ x ^ + ...). Besides the usual 
corrections to scaling due to irrelevant variables and the nonlinearity of scaling 
variables, this new type of corrections appears in the crossover scaling functions 
between two critical points since the matrix of anomalous dimensions cannot be 
diagonalized simultaneously in the ultraviolet and the infrared. The mixing is, of 
course, restricted by the symmetries that are preserved under the crossover, such 
as spin-reversal symmetry or self-duality. 

A multiple crossover involving the fixed points M p , M p _i, . . . , M 3 is character- 
ized by the appearance of p — 1 different length scales £ p +i, p , £ P)P -i, • • • ,£3,2 = £• 
At distances £ p '+i, p ',$r ^Cp',p r -i (corresponding to RG times v+i,p'~ t ~V,p'-i)' the 
behavior of correlation functions is governed by the fixed point M p /; for r££, they 
decay exponentially. At any such fixed point, the the expansion of the local magneti- 
zation contains the most relevant scaling field 0(2,2) of dimension ay = 3/(2p'(p'+l)) 
(the flow of subleading operators is more subtle due to operator mixing). Hence, in 
order to extract the leading thermodynamic singularity, one may integrate in Eq. 
(H) over the piecewise constant function x(r) = x p > for r p / +lp / < r < r p / )P '_i; this 
suppresses in addition the familiar corrections to scaling of the form (1 + Br' m ) due 
to an irrelevant coupling of dimension y^. The singular part of the susceptibility 
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then consists of a sum of terms 



■/£„/ 1 _/ \ Cn'+1.n' / 



2-2x p _i / t \ 2-2cc / 



Since x p > < 1 for all p', the leading singularity comes from the integration region 
£4,3 < r < the other terms are a third source of corrections to scaling. 

What makes this equation useful is the fact that any two subsequent crossover 
length scales have the same ratio 

= e 6 * . (8) 

This feature of the exact solution is consistent with the Lagrangian formulation 
(0), as can be inferred from renormalization group and scaling arguments (for details 
see U). Under the flow between M p and M p _i, the irrelevant running coupling u p 
and the relevant coupling m p _i mix, which can be expressed as the dimensionless 
equation 

ViV-T 1 ^" 1 = W p/yp ■ ( 9 ) 

Expressing these dimensionful couplings in terms of the crossover length scales 

V W-i ~ V-i ( 1U ) 

and using the exponent relation = A/p = —y p then shows that the ratio 
(£p+i, P /£p,p-i)/{£p, P -i/£p-i,p-2) is asymptotically constant. 
From (0), (H) and ([T0|) , one obtains 

7 = 1(2 - 2x p + 2 - 2x p _ x + . . . + 2 - 2x 3 ) = (P ~ 2) f P+1) . (11) 

The local energy density contains at any fixed point the most relevant field that is 
even under spin reversal, i.e. 0(3,3) of dimension ay = A/(p'(p' + 1)) for p' > 3 and 
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0(1,3) of dimension £3 = 1 for p = 3. Hence for the specific heat, one obtains in a 
similar way 

a p = — (2 - 2x p + 2 - 2x p _i + • • • + 2 - 2£ 3 ) = p(p ~ 3) ; (12) 
V P 2 

the additional logarithmic singularity in Eq. (||) comes from the integration region 
£4,3 < r < which is governed by the Ising fixed point M3. Zamolodchikov's S- 
matrix has also been related to multiple crossovers in the D-series of minimal models 
0, based on a thermodynamic Bethe ansatz analysis for antiperiodic boundary 
conditions ||. A £)-series staircase model has the same Lagrangian description (g) 
and the same flow of scaling fields that are even under spin reversal, but Eq. flTT|) 
is replaced by a similar sum over scaling dimensions of odd .D-series fields. 

From a renormalization group point of view, two (related) aspects of this system 
are remarkable. 

(a) The irrelevant parameter t p in (Q) affects the leading exponents associated to 
the relevant perturbation t p . This effect, encountered also in regime IV of the 
integrable lattice models of Andrews, Baxter and Forrester 0, is due to mixing 
under the renormalization group 0. 

(b) Several fixed points democratically share the responsibility for the exponents 



(|TT|) and (12|). Generically, the ratio of different crossover "times" T p / p /_ 1 — r p / +ljP / 
would be singular as t p — > 0. That this is not the case in the staircase model is yet 
another of the amazing fine-tunings occurring in integrable systems. 
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Fig. 1. The ground state scaling function f(p, 9 ) for the staircase model M(9 ). 
All steps have the same logarithmic width 9q. 
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